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ABSTRACT 

Attitude control techniques foQ? the pointing and stabilization 
of very large, inherently flexible spacecraft systems are investigated. 

The attitude dynamics and control of a long, homogeneous flexible beam 
whose center of mass is assumed to follox7 a circular orbit is analyzed. 

Jh. this study, first order effects of gravity-gradient are included, whereas 
external perturbations and related orbital station keeping maneuvers are 
neglected. A mathenatical model which describes the system rotations 
and deflections within the orbital plane has been developed by treating 
' ttie, beam as a number of discretized mass particles connected by massless, 
elastic stauetural elements. The uncontrolled dynamics of this. sy&l;sm 
are simulated and, in additiqni^l^e effects of the control devices are 
considered. The concept of distribtited mMal- eo which provides a 

. means for controlling a system mode independently of all other modes, is 
examined. The effect of varying the number of modes in the model as wall 
as the nunber and location of the control devices are also considered. 
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. NOMENCIATUEE 

= location of mass nig^ (or from the origin point 
of the main body 

b = location of mass nij^ fccm the origin point of the 

main body 

El = bending stiffness of uniform beam 

F = control vector (includes torques and for>ees) 

f =; scalar actuator variables; alsoj feedback control 

gains 

f = control force vector acting on the discrete mass system 

f ^ = residual coupling coefficent 

= principal moments of inertia of the main body 
K = stiffness matrix of the system ... 

L = length of the beam (L = 2 %) 

M = mass matrix of the system 

Jf . = mass of the main body 

m = mass of the end mass 

m. = generalized mass for mode i; also refers to end masses 

mg “ mass of the interior mass located in between and m 2 

N ■ • = number of modes 

P - ■ number of actuators 

' = modal coordinates 

R = radius of the orbit 

r^ .= radius vector from center of earth to mass, m^ 

T ’ ” kinetic energy of the system; also, control transformation 

matrix which provides transformation from discrete actuator 
variables to distributed actuator variables (f = Tu) 

V. ■ ' 



X 

x,s 







n 



$ 



*^1’ '*’2 





orbital energy of the system 

rotational enersr of the system. 

translation energy of the system 

torque acting on the main body 

elastic energy of the system 

independent distributed actuator variables 

potential energy due to gravity forces 

deflections of ’the end nesses Cm^ = m2 = m) 

deflection of Tiessessin^ and nij^, respectively 

state vector of the system (discrete coordinates) 

local vertical coordinates in the orbit plane (z - along 
the local vertical with the orig:m at the sjVstem center 
of mass) 

pitch angle of the main body 
eigen values (modal frequencies) 
orbital angular velocity 

eigen vectors (mode shapes) 

coordinate transformation matrix vdiich provides trans- 
formation from discrete coordinates to modal coordinates- 

relative angular motions of bo-di end masses relafive 
to the undeflected configuration of the beam 

coordinate system with origin at tlie interior point of the 
beam, (parallel to the x, z axes) 


vi 


KEST OF ILLUSmTIONS 


Fig. 2.1 System configuiration. with jnain lody and each beam 
nodelled an end mss 

Fig. 2.2 System configuration with min body nnd each beam 
modelled by two msses 

Fig. 2.3 Effect of gravitational forces on the system with 
main bocfy and each beam modelled by an end mass . . . 


Fig. 3.1 
Fig. 3.2 

Fig. 3.3 

Fig. 3. It 

Fig. 4.1 : 
Fig. 5.1 : 

Fig. 5.2 

Fig. 5,3 

Fig. 5,4 

Fig. 5. S 


Three-mass system configuration 


Possible approximtion of a firee-free beam by 
three discrete particles (first two modes Shown) . . 

Variation of and ui 2 with nig (local vertical 
system) 

Variation of co-, and with m^ (local horizontal 
system) — 

Ifode control inplemeirtation ..................... 

Uncontrolled motion of the three-mass system 

Cfi=0) 

Dynamic response of the system with the number of 
actuators equal to the number of modes (f^ ~ 1) .. 


Time history of the control forces with the number 
of actuators equal to the number of modes (f = 1) 

Dynamic respoi^e of the system with the number of 
actuators equal to the number of modes (f^ = 0.1). 

Time response of the system with the number of 
modes - 1) 



1. INTROrUCrXON 

In i±ie recently ccaipleted "Outlook for Space Study"^ a nuinbeip 
of proposed new space missions were shown to require large-scale, 
light-wei^t space structures. Three representative proposed future 
missions utilizing such ^sterns are: 

Cl) ocean data systems involving a 100 m. wide structure 
for the purpose of collecting data on the state of the 
oceans, pollution (both in the atmosphere and ■fee oceans) , 
and salini’ty. 

(2) electronic mail systoiis requiring a 50 m. diameto? antenna- 
receiver system to be placed dn a synchronous equatorial 
orbit. 

(3) a space-based solar power collector system, also operating 
in a synchronous equatorial orbit where ■fee arrays for 
collecting fee incident solar en^gy would have fee dimensions 
on fee order of kilometers. 

It is evident that a complete new technology must be considered and 
developed so that these s'truc^tures can be delivered into orbit Cusing 
fee Shuttle Transportation System), deployed, and feen fliUy asseEribled 
in a space environnent. Because of feeir inherent size, fee tes'ting 
of such systens in a ground environment is not practical, ^fadeling 
techniques and scaling algorithns must be developed so *feat fee performance 
of -these sys'fcans- can be accura-trely predicted prior to launch and assembly. 
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Tn many circumstrances, it will be necessary to control the shapes 
of the antenna or collector surfaces to within centimeters or even 
nuJlimeters by using a variety of sensor-actuator systems to be positioned 
throughout the flexible members. 

To address some of these unprecedented problems a special Industry 

2 

Workshop on large Space Structures was held at NASA Langley Sesearch 

Center in Feb. 1976. Among the principal conclusions ^vas that technology 

development is most critically needed in; the definition of large space 

structural configurations; iiproved modeling and scaling techniques; and 

the interaction between, the control system and structural dynamic responses 

It was also stated that some of these development techniques would utilize 

irproved conpiterized analytical models. Concern was also expressed over 

possible resonant interactions between some of the flexible structural 

modes (with frequencies greatly reduced as ccanpared with more conventional 

structures) and the frequencies associated with the attitude control 

systems. . Analytic areas requiring model development also include the 

precision detecmination of gravity-gradient farces and raomente acting on 

2 

such large structures. 

The AIAA Symposium oii Dynamics and Control of large Flexible Space- 
craft held at Virginia. Pol 3 rtechnic Institute and State Universily, 

June 1977 provided a review Of the state of -the art in this area. It 
reflected the tremendous; strides made in modeling and analysis of space- 
craft in tbe last two decades. Yet, the next two decades are likely to 
place a severe strain on the current ; i^ate of the art if some of the 
advanced concepts, such as solar po?/er satellites or space colonies, are 
to be converted into re^Lity. 

2 
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In this r'< 2 gard, this synroosiuiii raised ns inany questions as it answered 

and, in the process, it pointed to areas of future researcli. 

It was stated in Ref . 2 (Page 12) that to assure the availability 

of the large space structures tecihnology to support future missions 

a comprehensive research and development program must be defined. 

2 

Three of the objectives of such a program are: 

(1) development of active surface control techniques and 
systems vdiich can measure and correct surface deformations 
to within millimeters of accuracy . 

(2) development of attitude control techniques for the 
pointing and stabili 2 stion of vsry large, inherently 
flexible space structures. 

(3) development of analysis and simulation techniques which 

can extend subscale ground test experience to hi^-confidence 
predictions of full-scale performance in the space environment. 

The present study represents a preliminary contribution to tfcie accomplishment 
of the second stated objective. 

In this report. Chapter 2 deals mth the development of a nathematical 
model which describes the rotation and deflection of a long, thin, flexible 
beam in the orbital plane. It is assumed that the beam is represented 
by a nnin body and a series of discretized particles and the model pre- 

♦ U ........ , , 

vously developed by Meirovitch and Nelson is used as a starting pomt. 

This model is extended to include first order effects of gravity-gradient 
torques. 



The discrete models are developed \mder idle assiunptions , first , that 
the beam is represented by a massless rod connected by two end masses, 
and, secondly, that the beam is represented l:y four discrete paiTticles, 
two at the ends and two located at points taken half-way between the 
center of the undeflected beam and the end points . The equations are 
developed and linearized for the cases of small amplitude rotations and 
deflections. These equations are then cast in the proper form for the 
application of the mode control concept (Chapter 4) as eiqjounded ih a 
recent contract report by Rockwell Ibitemational. ^ 

A more general case of a flexible beam system using three discretized 
masses without a central body is considered in Chapter 3. The nonlinear 
equations of motion are developed using the lagrangian formulation. These 
equations are then linearized assuming snell amplitude defornations about 
two equilibrium positions : (1) alignment along the local vertical and 

(2) alignment along the local horizontal or orbit, tangent. Stability 
conditions for system motion taken about these equilibrium positions 
are obtained. 

The mode control concept^ for controlling a spacecraft by inde- 
pendently controlling motions of the rigid body and the vibrata®’^^ modes is 
prejsented in Chapter 4, Development of the mode control concept is based 
on two coordinate transformations . 

In Chapter 5 , the mode control concept described in Chapter 4 is 
applied to the linearized equations developed in Chapter 3 for the three- 
nHss system. The uncontrolled dynamics, as well as the dynamics of the 


systQtt wilii two actuators (equal to the raraiber of vilirational modes) , 
or one actuator, are simulated numerically using closed form expressions 


i 

The concluding comments from this study and suggestions for future 
H work are summarized in Chapter 6. 

I 

I 

I ■ ■ ■ ■ 
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2. MAIN BODY WEEB TOO 


2.1 Equations of Motion-Torque Ecee System 

2 ,1. 1 Each beam modelled by an end mass 

The configuration of the system consid^ed is shown in Fig. 2.1. 

The system canter of mass is assumed to move in a circular orbit and 
the system is constrained to move in the orbital plane. The system 
consists of a. main body 0%) with two beam? attached to ^ the main body. 
Each beam, modelled by an end mass , is attached to. the main body 
as shorn in Fig. 2.1. / The equations of motion for the system will be 
drived by use of the lagrangian formulation. One generalized coordinate 
of the system will be the an^e 0 xdiich represents the rigid body pitch 
angle . The other generalized coordinates, will be the deflections v^ 
and V 2 of the end mas.sesi m^^ and m 2 , respectivelyi -r^ativC to the 
undeformed Cbut rotated) body z axis. 

The total Idnetic energy Of the system, in terms of the rotational 
and translational energies, is 

T = T^ + T^ -f- const, 'due to circular orbital motion (2.X) 

The rotational kinetic energy of the main body is - 

'i- ' (2.2) 


1 *2 


6,7 


and the translational energy, due to the end masses, can be written as, 

2 ^ ' ■ ■■ ■: 

^ . .. (2.3) 

2W 


T. = f 2 CV. (SV. • £V.) 

I. 2 X X i oSf X X 


where 


^i = ^iW^ 


(2 id 4-) 


M >r t 2m 1 m, - nu = m>: :i 


(2,5) 




The coordinates of the two nesses and nig with respect to point 
0 (Tig. 2.1) are given in teznis of the position vectors, 


r^ “ ^ a- (2.6) 

A 

Fj = - Vgj - a k (2.7) 

After substitution of Eqs. (2.2) - (2.7) into Eq, (2.1), the resulting 
expression for the kinetic energy is 

T = 1 Ij' 8 ^ 


where 


r* ^ ^^1 " ^2^^ ^ ^'^'l " ^2^^ ^ ^ const. 


Ig + 2 na^ and m* = m^/M 


( 2 . 8 ) 


Initially the effect of gravity-gradient forces will be neglected so that 
the potential energy consists of the elastic energy, plus a constant due 
to the assumed circtOar orbital motion of Hie mass center. 


V = U t const. = ^ (v^^ + + const. 


(2.3) 


where k is Hie stiffness of each beam (k = 3EI/a^). 

The Lagrangian equations of motion for the system have the form 

(d/dt) (3L/3q^) - (8L/3qj^) = 0 ■ (2.10) 

where q^ assumes Hie values; 9, v^, and V 2 * The langrangian has 
the form 

. L X T - V = T.“U -J- const.: . . 

where T is the kinetic energy and U is the elastic potential ener^ of 
the. system. 


Gonsider±ng only tiie linear teGnns, the equations of motion can he 
represented in matrix form as . : 




C2.U.) 


If the coirtrol forces Cincludang the control torques also) acting on the 
system are represented hyj F, then Eq. (2.11) is reduced to the standard 
form . ■ 


M X t K X = F 


( 2 . 12 ) 


For the special case, when the deflections of the end masses are assumed 
to be antisymmetric (v^ - = v),Eqs. (2.11) and (2.12) reduce to the 

following equation: 




I r 

I ' 
-*-2 

2ma 

1 Q 

2ma 

2m 

TV ^ 


J 




2.1.2 Each beam modelled by two iiHsses 

The mathematical model .for ..the system considered now consists of 

a main body and two- flexible beams which are attached to the maiix body. 

Each beam is modelled by' two masses as shown in Fig. 2.2. For simplici'cy, 

it is assumed also that the elastic motion of the beams Is antis5?mmetriG, 

:■ that is V = V ' and v, = v-u' * 
a a b b 


WitJi 13163 e restrictions, the Jdnetic enef’gy expression takes the form 


1 *9 • • 9 *9 9 9 * m fp 

T = j Xg 9 + m [(v^ + 9a) + ) + CVj^+9b)^] 

+ const. 

The elastic strain, energy is 

" = %'^a^ * * ^22'^b^ 

Fca? small displacements, the equations of motion can be written as 


(2.14) 


(2.15) 


I 2 " 2ma 2irib 


2ma 2m 
2mb 0 


0 

2m 


L, 

where 


• • ^ 


v_ ( 


Vi. 


J 


1 


'I 


0 


21^ 


2k. 


12 


2ki2 2k22 


V- 


Vt. 


v= 


V, 


b J 


I2" = I2 ^ 2 m (a^ + hh 


(2.16) 


Eq* C2.16) is in the form of Eq, (2.12) ^ere X = (Q v, )^ and F = 

a. X) 

The stiffeess natrix [s] for the beam modelled by n nasses can be 
obtained from the flexibility matrix [a] using the relation [s] = [a]""^. 
■Hie displacements .Cv^) in terras of the flexibility influence coefficients 
(cij) and the forces (f^) xhich are associated with the displacements can 
be written as 


'^l ~ ®‘n ^1 “12^2 “in 

» t t 

= ‘il 


In mtrix: form, we can write 
(v] = [a] [f] 


(2.17) 


For iihe case ^^ere each beam is modelled by ti-jo masses, the stiffness 
matrix becomes (b = a/ 2) 


[s] = 


% 

k ^2 

^2 

^22 


El 

~:r^ 

7 a 


2 

-5 


2.2 Equations of Motion - Gr?avitatioml Effects 
2.2.1 Each beam modelled by an end mass 


-5 

16 


C2.18) 


The configuration of the system with main body and two flexible 
beams attached to the main body is shoxei in Fig. 2.3. Here, each beam 
is modelled by an end mass and the elastic motion is assvnned to result 
only from antisymmetric bending Cv^ - Vg = v). Ihe kinetic ener^ of 
the system can be written ; 


1 9*9 *9 9*9 •• 

T = ^ (I 2 + 2raa'^) 9 t m [ v + v'^9^ + 29 av ] • 

^ due to circular orbital motion C2.19) 

The elastic strain energy is reporesented by 

U = C2.20) 


The potential energy due to gravity forces is eiqiressed as 

V = v„ - im C-j^^ + i ) ; *ere u = 

■ q 

where the potential due to the main body is 

V - - | “ 0 ^ '^1 - 2 « 


( 2 . 21 ) 


and the radial distances of m^^ and m 2 the center of the earth 


are ' 

W 2 ) ~ ■’* " 2 3^0 a cos {tf - (9 + } ]^ 

\r^\ - [^ 0 ^ t ~ 2 r^ a cos C9 -f- 


10 



3 3 

use of the isLnomial theorem and neglecting terras of order ( ‘/rg) 
and. hi^er, one can dev^op for the gravitational potential, Vj the 
expression 

(Ii-l3) cos2 9-y^ 

z nir^ 

- W, 


{1 + 3 cos 2 (9 + 


a-j (2.22) 

lagrange's equations of motion can be obtained and then linearized as 


Ij' 2ma 


2raa 2m 



2 

6Wq raa 


e 2 

otJijQ na 
6wQ^m+2]^ 


9 


V 


'9 




(2.23) 


where 


- It + 2ma'' 


2.2.2 Each beam modelled by two masses 

The equations of motion, with the assumptions stated for 
the case where each beam is modelled by two masses (m^ m^^ = m), can 

be developed in a manner similar to that used in Section 2.!2*1. Here 
only the final matrix form of the linearized equations is presented: 


Ig 2na 
2ma 2m 
2mb 0 


2mb 

0 

2m 


9 


v. 


V 


■ - I3) 

Bwg^ma 

ewQ^rob 




CD 

C ^ 

Dojg ma 

2 

SWq m+2kT T 

^^2 

/ 

I'^a , 


\\ 

I 

CD 

E 

0 

^2 

•2' ■ 

SUq iii+2k22 


■V'-K 

J 


1 K 

V 

\ D J 


where 


I, " I, + 2m (a"+ b^) 

1 J. 


2^ ^2, 


(2.2L^) 


n 


3. IHREE-MASS SYSTEM 


3.1 Eqxmlioiis of Moldon-Local Vertical 

The long beam modelled by three masses as illustrated in Fig. 3-1 
consists of two end masses (m^ = m 2 = m) and an interior (point) mass , 
rag. The system center of mass is assumed to move in a circular orbit 
and the system is constrained to move in the orbital plane. The' equations 
of motion arc derived ' using the lagrangian formulation. The generalised 
coordinates and ^2 ^:^pa^e3ent the relative ■’ti^ilar motions of both 
end msses relative to the undeflected orientation of the beam. 

3.1.1 Expression for Idnetic energy 


The total kinetic energy, of the system can be written as 

T = + T^ + T^ (3-1) 

The rotational and orbital ^ergies are 

= Oj T^ - ^ (3.2) ' 

■ B 7 ■ 

and the translational energy is obtained as ’ 

1 ■ ' ' 

T . = i m 2 (y. ‘ V.) - 5^ (2 V. * 2 V.) . : (3.3) 

L a i i “ 

where M = 2m + mg , ojg orbital angular velocity and R = orbital radius 


The velocities V^and V 2 can be obtained from Eq. (2.4) as (Fig. 3-1) 


^1 ~ ^y ^ ^ ^ 


^ ^ tUy j X (^ 2 ^ sin^Jji- + ^ 2 _ cos k ) 


(3.4) 


12 


Vj = j X C- Ij i - Sjk) 


= 03y j X c- Jig ~ ^2 


(3.5) 


where 


CO 


W = <i>ii V = “o ^ '*’2 


(3.6) 


1 " ^ •>'2 

^ substitution of Eqs. (3.2) - (3.6) in Eq. (3.1), we find that 


T = I'M* [{(coq t ^^) sin(J)^ + (co^ + ^g) Ag sin4>g}^ 




+ iHq {(tOg t 1^)^ A^^ t (cOq + ^g)^ ^2^^^ 


-f- M R^cjp^/2^ 


(3.7) 


_ 2 ,^ 


where M* = m'^/M, the systm reduced mass, and 5 ^= 111 ^/ 111 . It should be noted 
that the last term in Eq. (3.' 7) is a constant for the case of a circular 
orbit and does not affect the equations of motion. 


3.1.2 Expression for potential energy 


The potential energy V of the system in an inverse square force 
field is given by . 


V = -w{|fT +-S- 


^ll 1^2i kol 


> 


(3.8) 


where @1 is the earth? s gravitational .constant , ,a3^ \v*\ is the distance 


of from the center of the^herical ear^. Tien |r^| may be expressed, 
in terms of R = |Rj and the coordinates of in the local vertical system 


with origin yat the center of the beam, as 

2 


2 


Pll = :<R + ?i + Com> * ^ 1 




(3.9) 


IF,I = [CR- 5, + o -i- ( -So + 5™.) ] 


C3.10) 


l>^ol = I'K * ^ 


(3,11) 


where the coordinates of the system center of nass measured with respect 


to point j m_, on the beam are given by 


«om = - " <«1 - ’ f'l- 


(3. 12) 


The expression for V, after omitting the terms of order 1 / R and hi^er 
in the bincanial expansion of l/|r^l, becomes 

V = M* *’* ^2 sin$^ + ^2 sinil)^)^} 

+ mjj {(^2^^ cos^(j>^ + S,^ cos“(j)2) - ^ (S^ sin^({)^ t sin^4'2^^] 
3.1.3 Expression for elastic energy 


The elastic energy’'*^ is obtained by assuming that the end masses 
move as cantilevers with respect to the reference point, mQ, on the beam 
(Fig. ; 3.2). Thus the elastic energy in terms of elastic defcfrmations 
and expressed as 

;:U= |[k3_ 5^2 t kg 52^] 


S.^ sin-iji^ + kg six? 


(3.14) 


For the later appxodmate modelling of the elastic deflections of a free- 
free beam it will be convenient to assume that = Ag “ ^ such that mg 
and the system center of mass will be coincident xdien - 4>2 ~ 0 (see Fig, 
3.2);.- -■ 


3.1.4 Lagrange's equations of motion 

^ The general equations of motipn are- developed using the lagrangian 
formulation for the variables. i = 1^2. The Lagrange's equations 


are 


d 3T 3T . 3V . 3U _ 
dt ^ " 3(j)^ 34 k 94)^. “ 


(3.15) 


^ 1 . '-L 

The equations of motion are dbtaijned, using Eqs.> (3^7), (3.13) ‘and (3.14) . 
in Eq. (3.15) , as 

«« . *♦ 

M* [{4>]_^2. \ cos^^ t 4>2 ^2 ^^*^2 

t (tijQH2) ^2 ^ sin$^ 

+ {(Wq+4|j^) a^sin4>^+ (tiiQ+lg) ^2 cos4)-,4'^ 

+ { 4 ^ cos4]_ -(Mq+ 4)^) S,^ sin4)^ 4>2^ + 42 ^2 

~ (Wq+ 42^ ^2 *^*2 ^2^ cos4^ 

- {(Wg+4^) ■*■ ^2 ^1 

r " , 2,, 

. TT ^ 'h. *1 

-M* [{((Ugt 4 ]_) % sin 4 ^ + (Wg+42) i!,2 ®^'f'2^ \ cos 4 ^ 

- {(Wg+I^^) cos4^ + (tOg+ 42 ) gos 42> sin4^ Cwg+4^)] 
-i- ZoigSl* cos 42_ ^>2 '^°®'J'2^ sin4^ 

+ i sin4^ t 5-2 sin423‘ cosd^ 


+ sin4^ cos4-j^ = T^ 


q 


(3.16) 








^ ii 




r' 


M* [{4>2_ Jl-L sin(t)^ .+ ^1 ^1 ^ ^2 ^2 


+ ( tiiQ +^ 2 ^ ^2 ^ 2 ^ ^<2 ^^*^2 

+ {(Wg+4i^) sincj)^ + (“0+^2^ ^2 ®^^2^ ^2 ‘^°S'i>2 ^2 


+ {(j)^ cosiji^ - C(j)q+|^) 5,^ sintj)^ + <J)2 5*2 costjj^ 

- (o)q+|2) ^2 ^2^ ^2 °°®^2 


- {( Wq +^^) Jt^ cos ( j )^ + (Wq+|2 ^ ^2 '^°^‘f’2^^2 ^2 


* T ‘*2 


- M* [{(Uq+I^) sdxi4i^+(wQ+l2V ^2 ®^^2^ ^2 ‘^°^'f’2 


- {(Uq+I^:) j!,^ cosd)^ + (Wq+$2^ ^2 ^2 

+ 2 uip^* [{a^ coS(|i^ + 2.2 Goscljg} 5*2 ®S^'1’2 
+ ^ sin4i^ + fi-2 siruf>2} ^2 


®0 . 3. ■ 2 ... 


+ — C-^) sin(|)^ cosA^] 


+ IC2 sfa*2 cos*2 = T^ 


C3.1'; 


3.1.5 Linearised equations of motion 


Eqs. (3.16) and (3.17) are linearized about the local vertical for 
small amplitude motions such tihat sin(j)_* « and cos*. ~ 1. Also as 
indicated in section 3.1.3, we now assijme, in order to model a free-free 
beam, A, = = I and Ic, = 3<u = k. 




B7 letting = v^/il and ~ ^ 2 ^ linearized equations of motion 

can be represented as 


M* (1+ mg) M* 

M* M*(l+inQ) 









M*C2-HHq) 


+ k 



3Wp^M*C2+mQ) 


+ k 




The control forces that are assumed to act on and m 2 are represented by 
F F 

and V 2 J respectively. For the special case of a two-nass system 
VTith jiIq = 0, Eq. (3.18) reduces to ■ S 



The characteristic equation for this two-mass system is obtained as 
s^ + (3 wq^ t k/m ) = 0 

This eqiiation shows that . the system is forced to oscillate in an anti- , 
symmetric mode if the ' flexible beam is modelled only by the two end masses . 

■' 'I " 

3,1,6 Stability ■ an^ysis ,f 


The characteristic equation for the three-^ss system from Eq. (3,18) is 


17 








( 3 . 


CM*) (2mQ + mg“)s 
+ 2M* (1+inQ) {3 uq^M* (24mg) + k}s\ 

+ {3u)g^M* (2+1^) + k }^ = 0 


The magnitudes of the two natural jErequencies , and oigj as a fimction 
of the central mass, ra^j with M = 2m + = constant, as obtained by 

solving Eq. (3.19), is given in Table 3.1. 


i3(iiQ^+(3k/M)]*^ [3{3uig^+(3k/M) }]^ 


[3cjQ^+(2k/M)]*^ I [2{3t0g^+(2k/M)}]^ 


[3o3g^+(5k/3M) L|{3t0g^+( 5k/3M) }]' 


[3/jjg^+(3k/2M)]^ [| {3u)g^+(3k/2M)}]' 


Tfeble 3.1 Variation of and with nig 

; The following system parameters and initial conditions are assumed 
for numerical study: L = 2A = 100m; orbital altitude - 463 km (250 nautical 

miles); El = 7.707197x10^ N-m^; k = 3EI/Jl^; M = 1000 J<g; v^(0) - O.Obn and 
V2(0) = v^^CO) = V2(0) = 0. Structural pasrameters are taken for a beam 
made of wroiight aluminum (2014 TB) and c yl i ndr ical in shape. The outer 
diameter of the beam is 50ram and the thickness is 5mm. The variation of the 
two natural frequencies, tod with :Hg is shown in Fig. 3.2. Also shown 
in the figure are the variations of and CO 2 with nig when k = 0 , indicated 
3^ £Uj_g and ci) 2 q. 
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We observe that remains constant at ^/^ Wg (which is the frequency.. 
at which a rigid, dunibell satellite will oscillate) , but CO 2 Q decreases 
with an increase of iHg. For the general case, increases with mg, 
but W 2 decreases up to a certain v^ue of but then increases with 
ar^ increase of iHg. 

3.2 Equations of Motion - Local Horizontal ' : . 

3.2.1 linearized equations of motion 

In order to develop the small amplitude equations of notion for the 
case where the beam is liominally al igied along the local hcrizontal , i^e 
begin by replacing ^ Eqs. (3.16) and (3.17) by (tt/ 2 + (j)^). After 
appropriate simplification and subsequent linearization the equations 
of ^tion may be represented in mtrix form as : 




- {3Wg V (1+nig)- k} 

- awg^i* 




’’ L-' 


i 

. = 1 i 


- 3o)g^M* 

- {3u)gV (1-Hng)- k} 

l.V^! 

CM 


(3.20) 

For the two-mass system with nig - 0 Eq. (3.20) reduces to 
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The charac±:e3?'istic equation for the two-irsass system is 
.+ (k/mV - = 0 

The system behaves like a duiiibell for k > 3 (i)q m and is unstable for 

„ 2 \ ^ ^ ■ -r- 

k ^ 3(0p m. - V . ' 

3.2.2 Stability analysis. , - 

The characteristic equation for the three-mass system can be - 
developed from Eq. (3.20) as 

+ [BCM*)^(Dq^ - 2M* Cl^){3Wg^M* (Irh^ 

. + [{3WqV a+mQ> - k}^ - 3(.1 q^ = 0; ; . ^ 

The variation of and W 2 with^miQ is shown in Table 3.2 with M :- 
2m + mu. = constant. 


[C3k/M> - 3ti)Q^]^ 
I(2k/M) - 3 Wq^]^ 


[3{(3k/M) - Wn“}] 


[2{(2k/M)- liSWp^}]^ 


[C5k/3M) - 3Uf, ]' 


[| {(5k/3M) - 


[C3k/2M) - 3 u)q- ]' 


r| {(3k/2!0 - 


■Jable 3.2. Variation of and wi^ with mp. 


The variation of the two ratural frequencies and W 2 jibr the three-mass 
system with tUq is shown in Fig. 3.3. For the assumed system numerical ^ 
parameteus, both 01 . and tOo decrease with an increase of m^. ' 








4. MODE CONTROL CONCEPT 

A cmncept is pcesented in a recent Rodcffell Intennaiional Report 
(Ref. 5) for controlling a flexible space structure by independently 
controlling iKotions of the. structures ’s rigid body and vibrational 
modes . The mode control concept leads to criteria for locating 
actaators and algorithms for their combined use. Development of the 
mode control concept is based on two coordinate transformations. Dis- 
crete displacements and discrete forces are transformed to modal coordi- 
nates and distributed actuator variables, thereby completely uncoupling 
the system equations. In this chapter, the mode control. concept is 
outlined briefly for our specific application but more complete details 
are available in Ref. 5, 

4.1 I)ynamic Equations 

A flexible space structure is modelled as many rigid bodies inter- 
connected by massless, elastic structural elements. Then small amplitude 
motions are descm»ibed by the linear differential equation 

M X + K X = f^*^^ ' (4.1) 

^diere X is an N x 1 vector of discrete coordinates measuring the angular 
and translational displacements of each body relative to its inertially 
fixed rest posi-tdcn, M is the real N x N. mass matrix, K is a real, 
symmetric N x N stiffness matrix, and f^'^^ is 'an N x 1 vector of the 
control forces. 
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4 . 2 Eigen-Analysis 


Hie eigenvalues (modal frequencies)} } - — » o^} and associated 
Nx 1 eigenvectors (mode shapes)} ~ tj)^^^}are obtained from 

the homogenous part of Eq. (4.1). The orthogonality properties of 
the modes provide that 




0 i 7^ j 


1=3 




i j 


m^to. i - j 
11 


where is ref erred to as the generalised mass for mode i. Next, a 
coordinate transfooonation is performed utilizing- an N x N matrix, 
constructed feom the N eigenvectors 


A transfoimiation is introduced by assuming that a given displacement 
profile may be ejqressed in terras of a series of *Qie shape fimctions 
(mode shapes) multiplied by time-dependent weighting factors (modal 
coordinates i.e. X q^ + q^ + *-* + q^, or 

X = $ q "" ■■ ■■ ■ ■ ( 

- T 

Substitution of Eq. (4.4) into Eq. (4.1) and premultiplication by ^ 


produces 


M $ q + K $ q = ^ f^*^^ 


I 



Finally, the transformed equation of motion is obtained with the aid of 
the orthogonality conditions, Eqs. C4.2) and C4.3), as 


^ ^ *4 q = r <1 ■ 


C4.6) 



The advantages of the modal fonoolation are that the left-hand side 
of each scalar equation associated with Eq. C4.6) is uncoupled in q^ 
and criteria can be easily developed for . truncation purposes ^ Generally, 
f^°^ has low frequency components. Thus a few low frequency modal 
equations can be. used to predict the response of a d^?namic model, which., 
has been modelled initially with a larger number of discrete coordinates..; 

4.3 Actuation ^ ^ ^ • 

4.3.1 Nuniba? of actuators equal to the. number of modes 

Ihe discrete control forces, f^*^\ are transformed into 'generalized 
control forces, by 

fCg) _ p, ^/m^ C4.7) 

Independent actuation of all of -the N modal equations can be achieved if 
N actuators are used in such a way so as to ptroduce a generalized force 
in any given mode without, forcing the other modes . ■ Then we obtain 
f^^^ u if we let 

f*'^^ = M ? u V ■ C4.8) 

where u^. represent the independent actuator variables. After substitution 
of : Eq. C4. 8) into Eq. (4,6) one arrives at 

■ • : " . 2 

q 4 q = ,u • (4.3) 
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The niode control implementation for this case, where the ntmiber of inde- 
pendent actuators is equal to the numher of modes of the system, is shown 
in Fig. 4.1. 


4.3.2 IJumber of actuators less than the number of modes 

In most instances , independent control of all the N mode's ijf a 
dynamic model using N actii^prs is impractical since N is usUii';ly. ^ 
very large number. Additionally, control is often required of only a few 
of the lower frequency modes. For these reasons', it is n^essary to modify 
the foregoing procedure to establish a means of independently controlling 
only a selected small number of the systeiaVs modal coordinates, r 

Hare an expample is illustrated taking the number of actuators CP) 
as 2 and the number of modes CN) of tiie system as 3 such:. that P < N. 

[The general theory is available in Ref. S.] The control force for 
this case can be written as 

i: 


' f, 


rCc) _ 


•s 


r Cc) 

■2 

rCc) 


1 

0 


0 

1 


u 


■2 J 


F f 


(4.10) 


where and f^ represent the scalar actuator variables*; The generalised 
control for cev Eq. C4.7), is 


f “/m^ j F f 


(4.11) 


A transformation of the form 


f = T u (4.12) 

is required to produce independent generalised control tbz^es. 


Here T -provides a iaransformation, from discrete actuator variables to 
distributed actuator variables such that u^ produces independent actuation 
of the i^ modal equation. Eq.s. C4.ll) and C4.12) are combined to yield 

^"^rTu (4,13) 


Eq, (4.13) can be written, with N = 3 and P = 2, as 


r f Cg) N 
1 

'r- 

^/m^ 0 0 


‘*’11 ‘**21 ‘*’31 

Cg) 

1 ^ 

= 

0 ^/m2 0 


‘**12 ‘**22 ‘**32 

j 

i 

r 

0 

0 

CO 

\ 


***13 ‘*’23 ‘*’33 

** 


1 □ 

0 1 

0 0 


'll -12 


'21 ^22 


u. 


u, 


(4.14) 


Ejroanding Eq. (4.14), we obtain 


< f. 


Cg) ^ 


Cg) V 


^^3 


Cg) 


"11 * hi 


/^ll ^21 "^22^ 


m2 ^‘^12 '^ll ^ ‘*’22 "^ll^ 
^ *‘‘*’13 ^XL ^ ‘*’23 "^ 21 ^ 


^2 ^‘*’l?i "^12 '*’22 "^22^ 

^ ^‘*’13 "^12 ‘*’23 "22^ 


(4.15) 
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The traxisforsiation mtrix, T, is obtained by noting that for independent 


control of the first two modes,, we need 

f f 's’ 


(4.16) 


After a comparison, of Eqs. (4.15) ^d (4.16) for and f 2 ^®\ we 

can obtain T. Kiowing T, the residiWl coupling coefficients, and 

for the third mode can be obta^ed by .rewriting Eq. (4.15) using 
Eq. (4.16) as i, \ 


(4.17) 


V¥2> 


Thus, f ^ and f^2^are o]rtained frO^ the lasmrow of Eq. (4.15) after 
the deteniiinatian pf T. i A 

An interesting (and useful) resxilt bf this eJbniple is that, for 
the general case, only modes i=l to P are 'needed tcv determine T, and 
only and T are needed ip obtain the residual coupling elements for 
q^(i = P+l,..., N)i ?'^Thus, f;OV the case wherd'v^P<N, independent control . 
of modes i = 1 to P'|i5' possible and the respon^p of th|^ modes, N>P, 
depends on the residual coupling due to the P ac’Uiator^i. The mode control 


concept described in this chapter will be applied '|:o a jijong beam in space 
modelled by the three-^imass ^stem with the number of acfuators equal 
to or less than the number Of modes of the system* % | 


26 


5, MODAL CONTROL OF THE THREE-MASS SYSTEM 


5.1 Three-M^s System-Local Vertical 


The modal control of the three-mass system (Fig. 3.1) is now - 
considered. The linearised equations of motion [Chapter 3, Eq. (3^18)] 
are 


a 

h 


b 

a 




V. 


2i 


G 

0 


0 

c 






1^2 


Ji'F, 


Vr 


where 


a = M* Cl+mn) a b = M* 


G = 3 Wq'^ M* (2+h1q) + k 

The uncontrolled d 3 ?namics of the system is considered first and then the 
mode control concept outlined in Chapter 4 is applied to obtain the 
controlled system response. 


5.1.1 Uncontrolled motion 


Using the Laplace transform method with F = 0, Eq. (5,1) can be 
written as ' ^ 

P 2 - 

r s +ca 


v..(s) 




V2(s) 


cb 


sv^(0)tv^(0) 


(s^+ca)^- (cb)^ 


_ 2 — 
cb s +ca 


sv2(0)tV2(0) 


where 


c = c/(a^-b^) 



ISie solution to Eq. (5.2) is obtained, under* the assun^Jtion ihat 
(0) = V 2 (0) = Vg :C0) = 0, but (0) ^ 0, as 

y^(t) = i [cosjoj^t + cos cu 2 t] v^(0) 


(5.3) 


V 2 (t) = j [cos, w^t - cos £»i 2 ‘t] %^0) 


(5.4) 


where 


1^ .^*T^ 


0 )^= [c/(a+b)]- ; U )2 = [c/(a-b)]‘ 

The motion ox the center of the beam, iHq,, is obtained frcm Eq. (3.12) , 
and noting that = A sin = v^, 

- (mM) cos ti) 2 t v^' (0) ; =0 (5.5) 

The uncontrolled motion of the system with the assumed system para- 
meters and imtial conditions stated in Section 3.1.6 is obtained as 

v^(t) = 5 (cos 0.023635t + cos 0.040937t) mm (5.6) 

(5-7) 


VgCt) = 5 (cos 0.02363St - cos 0.040937t) mm 


^cm “ “ ^ 0.040937t mm 


(5,8) 


The time response of the beam end deflections, and v^, is shorn in 
/ ■ ^ 

Fig. 5.1. When hiIq = 0, the deflections and V 2 are the same in^cating- 

that the system.behaves like a dumbell satellite. The presence of ra^ 

produces the sScond frequency, co^, and the deflections due to are 

superinposed qh and v, as seen from Eqs. (5.3) and (5.4). 


28 


5.1.2 Number of actuators equal to the number of modes 

Considering the homogeneous part of Eq. (5.1), the eigen values are 

CS.9) 


2 2 

= c/Ca't-b); ^2 ” ~ c/(a-b) 


The coordinate transforaiiation matrix, is determined from ,4iie. eigen- 
vectors as 




1 

1 


1 

-1 


(5.10) 


The generalized mass matrix, from Eq. (4.2), is 


m.£ M 5> = 


2Ca+b) 0 

0 2Ca-b) 


C5.ll) 


The modal equation given by Eq. (4.9) for this case becomes 


‘' 1 1 


c 

( 

^-1 

( 


F -}. 




■ = 1 1 



0 Xj 

1 

^12 

V. -J 

1^2 } 


(5,12) 


We 


assume a proportional displacement and rate feedback for u^ and 


U 2 in the form; 


“1 = - % - ^2 % 


“2 ' ta °-2 % 


(5.13) 

(5.14) 


Ttie equation for q^ ;dLth control u^ is written as 
^1 ^2 % ° 


(5.15) 


The solution of Ea-t C5^15) as 


£« t -F 

2 ^ “2 1 

q^(t) = e [cos w^ 2 't + ~ — sin 


“12"-' ^1' 


CS.16 


where 


'^12 ” *’* = v^C0)/2 


(5.17 


similarly. 


where 


_ -u 

■ ^ 4 1 

q 2 (t) = e [cos — sin Ugi^t] q 2 ( 0 ) 

34 


“34 "" t^2 + fg - Cf.^^4)]^ ; q^vO) = v^CO) /2 


(5.18) 


(5.19) 


The discrete coordinates, and are related to the modal coordinates, 


q^ and q 2 , using Eqs. (4.4) and (5.10), 


'^1 = * V ''a = ii - <12 


(5.20) 


The control forc^ are obtained from. Eq.- (4-8) as 

T ^5? V,_ + (*1.^) 


+ (-2-i-f4) V, t (-2- h ) v„] 
3 3 ■ 




+ (^2-£4 ) V-, + ) v,l 

3 ^ 3 ^ 
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As a special case, when £2 = Xj^ and also, ” “ 34 * note that 
q^(t) = that the displacement V 2 Ct) = 0 and = 2 q^ for 

the assumed initial conditions. This condition is achieved if f ^ + 

- X 2 and f 2 = f^j^. Thus, by properly selecting the feedback control 
gain s, it may be possible to control a portion of the beam such i±at it 
will not be subject to any deflections. 

The uncontrolled motion of the system, obtained by setting f^ = 0, 
i = i, 2,3,4 in Eqs. (5.16) and (5.18) and in tum in Eqs, (5,20), is 
illijBtrated in Fig. 5.1, The dynamic response of the controlled system 
■wiih f^ = 1 is shown in Fig. 5.2. It is observed that 1he tip deflection 
amplitude, v^, is reduced to 0.01 urn from an initial„y^ue of 10 mm within 
12 secs. The initial control forces F,. and F„ are calculated to be 
-2.22N and—l.UN, respectively. The time history of the control forces 
is shown in Fig, 5.3. . The initial anplitude of the control forces can 
be reduced by reducing the feedback ^ins. Fig, 5.4 illustrates the 
time response of the system with f^ = 0 . 1 , a value which is 1/10 of the 
previously considered value for f^. In this case, the time required to 
reach a deflection amplitude of 0.01 mm from an initial value of 10 mm 
for the tip deflection v,, is about 120 secs. Thus a reduction in the 
control force increases the time constants of the syst^ proportionally, 
as Qcpected for a linear system. 

5,1.3 Number of actuators less than the number of modes 

Following Section 4.3.2, we obtain; the equations for q., and 02 as 

q-j_ + Xi q^ = u^ (5.23) 

q2 ^2 ^2 “ u^ , (5,24) 


31 


Eq. (5.23) can be controlled independently as in Section 5.1,2. Using 


= -f^ (q, - f 2 q^j Eq. (5.24) becomes 
*^2 ^2 ^^2 ^^1 *^1 ^2 

where 

g = (a+b)/(a-b) 

(Nolfe as long as mQ>0j a>bj see Eq. (5.1).) Using Laplace transform 
techniques, the response of the tiode q 2 due to the residual coupling of 
the actuator is 







For large t, is completely damped out [Eq. (5.16)] but 0,2 ("t) oscillates 

at the fbequency, m 2 C^q* (5.26)] . | ' 


The control forcej is obtained as 

= - m [Ft^Cv^ + V 2 ) + f 2 (v^ + v^)] (5.27) 

The time respon^^& of the system with the, ruiriber of actuators Idss than the 
number of modes is shown in Fig. 5. 5. The figure shows that the v^ response 
becomes oscillatory as time increases Jeven thou^ its initial value \®s 
zero. The maximum amplitude of the oscillation for v„ reaches ;-a value 
of 10 mm. It should be noted, that as time increases both v^ and V 2 

-t ;■ 

oscillate at the same frequency - i.e. that of the (second) uncontrolled 
mode- 


6. CONCLUDING COMMENTS 
6 . 1 Summary of •the Conclusions 

■fe a. consequence of the present analysis and numec'ical jresultSj 
the following conclusions can be : nade; 

1. The equations of motion are developed for a large space S'tructure 
system consisting of a main body mth two attached flexible beams, where 
each beam is modelled by an end nass or by two masses, CIt is assumed 
•that •the system center of mass .moves in, a .circular .orbit and the elastic 
displacement of idle beams are anti-symmetric,) It is seen that the 
s-tiffiiess matrices will contain contributions due tp both, the elas'tic 
strain energy coefficients and also the gravitational restoring effects. .. 

2 . The three-mass system is stable about the local vertical for : 
small amplitude motions but it behaves like a rigid dumbell satellite 

when "the central nmss is removed. It is . also found that the two-mass . 

1 * 

system is stable about the local horizontal system if;.the stiffness of 
•idle beam is greater than a certain value of •idle restoring effect pro- . ^ 
educed by "tdic gravitational forces on the end masses Oc > 3Wg m), 

3. Closed-form ahalytical solutions- are bbtained C'the local vertical 
system) for •idle tip displacements of the long beam which is modelled by 
the three masses, 

4. The UDdal con-fcroT concept is Very useful for independent controi 
of the system’s lower frequency nsodes . ‘ 


For 1±e tihreeriiass system where only . two mcdes are present, it is 
possihle to keep one portion of the beam without any displacement by 
properly selecting the feedback control gains for the. assumed initial 
conditions of the system. 

: '5. When the three-iiBss system contains only one actuator^ the 
modal control produces a COTvplete control of the first mode, but 1lie 
sedond mode isi not controlled. This indicates the disadvantage of 
the mode control concept when all the modes can not be completely , 
controlled., ^ 

,6.2 Recommendations for Future Wbrk 

As a result of this investigation the following topics are suggested 
for future 

1. In Chapter 3, the loiig beam is m,odelled by three discrete masses. 
For a better und.erstanding of the system response, the modelling of the 
beam by a greater number of discrete mass points should be considoced, 

2, The study could be extended to consider idle three dimensional 

motion of the beam by following the approach described in this report for 
• ■■' 1 : 
planar motion, ' - 

3, Optimal control theory could be used to obtain control laws 
by minimizing cartain cost functionals involving the state variables 
and elements of the control vector, 

4. The mode control concept could be applied to consider the control 
of a continuum model of the beam and the system response obtained from this 
model GOul,d be compared with the system response obtained here for the 
discretized model of the long beam. 

■ ^ 35 



5- , The deformatibn ox -the t)eam due to externally induced solar 

radiation effects should be considered both during steady-state and 

controlled responses. 
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Fig. 2.1. System configuration with main Isody and each 
team modelled by an end nass.; - . 





Fig. 2.2. System configuration with main body and each 
beam modelled by two masses . 





^ — actual Learn 



SECOND MODE 



Fig. 3.2. PossdLle approxiraation of a free-free beam 
by tiiree discrete particles (first two modes 






M =, 2m mp 

























5.4, Dyiiandx: response of 1ihe system V7x1ih 3^^ 
of actuators equal to the number: of modes Cfjj^ 
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uarors less than the number of ncdes Cf. 

^ T 




COMPUIER PROGR/yyl 


- -rOOB QCTAjjyj. 


1758/78 lOiSQjq 

IJOB [READ I# AT. 10:^19:551 SELLAPPAM 

IF0RT/A/8/E/P/S- F0RT.LS7L 

aiSTlNG ' ■ 

C, — FORTRAN STATEMENT 


c modal control 

EXTERNAL RGSOlfRGSO? 

dimension PARM(51,ycai^DX(/}3,Sl2ECa)fWORK(8,a3 

COMMON Fl,F 5 fF 3 ,F«,G 
CALL IN 0 UT( 5 / 5 ) 

CALL OPEN ( 1 , *SELLAPPAN" , 3 >IER 3 
.IFCIER.NE.l) STOP UNABLE TO OPEN FILE 
REaOC 5 , 91 ) TMAX>STEP,T 0 L 
- READ(g/ 9 l) X ■ . 

REA 0 C 2 / 9 n SIZE 
READ( 5 , 9 l) Fl,F 5 ,,F 3 ,Fa,G 

91 FORMATC 8F10.0 3 

PARM(t)s 0,0 

:PARM( 23 s.TMAX 

: P'ARMC3TgSTEP . . . ' ■ - . . . 

M=iJ 

WRITECS^Pa) TMAX,STEPrTOL 

92 FORMATC * ITMAXs %F8 . 5. 1 OX r * 5TEP= S F8 . a* 1 OX , ^ TOLs * f F8 . 6 3 

RKSCLCN,SIZE,DX,T 0 L,:PARM 3 
CALL RK 6 S (FARM, X,DX,N,IHLF,RGS 01 ,RGS 02 ,W 0 RK 3 

' ' : ■ . , . 7 ;: : ^ 

99 FORMATf *QIHLF=%I33 
GALL EXIT ■' 

END 


COMPILATION SUCCESSFUL OBJECT CODE IN FILE NAMED OOI.PR 


1 F 0 RT/A/ 8 /E/P/S F 0 RT.LS 7 L 
iLISTING 



C ^ FORTRAN STATE^EMT 

subroutine ROSOl ( TrVrDX) 

DIMENSION xca) , OX (4) 

COMMON F1,F2,F3,F4,G 

Dxm=x(gi ^ 

W0s0,001 US 
AK=0, 0005549 

AL1-380*^0*W0+AK ' . ■ 

~ DXC2 3=-F2*XC21-'(:AU+F1)*XCU 
DX(3)=XC43 
AL2"3 , 0*AL1 

DX C4 3 “^F4*X C4D » f AL2+F'3:) *X C 3 ) 'G*F2*X C 21 ^G*F1 *X C iV “ 

RETURN 

■ ;.ENb:-^.. 


COMPILATION SUCCESSFUL— OBJECT CODE IN FILE NAMED 002. P6 
JFORTVA/B/E/P/S F0RT.LS/L“ ' “ ____ ___ 


Ca . * . FORTRAN STATEMENT 

^ ^ ^ ^ RGS02CT>XrOXfIHLFVN,P) 

LOGICAL RKNXT 

DIMENSION X(4) ,DX(4) ,DUMMYC4I 

common F1,F3,F3/F4/G . ■ 

; CALL RGS01(T»X,OUMWY) 

Uls.Fl*X(l)-F2*^(C2) 

U2g-F3:»rXC3)"F4jfX(4) - 

AMsIOOO, 0/3.0 
FVi=AM*cUl+CU2/3.0J 3 

FV2sAM*(Ul^fU5/3.0)3 

Vl^XCn+X(33 
Va=XCl)-XC3D . 

V3^XC23»XC43 • ■ 

V4sXC23»X(4) 

VCMa*C2,0,/3,03*XC33 
IFCnNQT.RKNXT(IHLF33 GO TO 6 

wRITEC 5,13 T,XC1),XC33 »VlrV2,VCM,FVl/FV2,V3/V4 
1 FORMAT fix, F9,1,9F13.7} 

WRITE RINARYC13 TrFV2fFVl,VCM,Vg,Vl,X(33/XC13 

S CONTINUE 
RETURN 

END - 


COMPILATION SUCCESSFUL OBJECT CODE IN FILE NAMED 003.RB 


iRLDR TMP/3 SYS/E 001 002 003 DPOjSSP.LB I-0RT.L8 
ILOADED 

<DELETE/V SELLAPPAN ; 

DELETED SELLAPPAN 
iCREATE SELLAPPAN 

t fxec 

— — OEiGlNAL PAGE Ife 

Gr POOR QUALITY 
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